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ABSTRACT 

The growth of Jovian mass planets during migration in their protoplanetary disks is one of the 
most important problems that needs to be solved in light of observations of the small orbital radii 
of exosolar planets. Studies of the migration of planets in standard gas disk models routinely show 
that the migration speeds are too high to form Jovian planets, and that such migrating planetary 
cores generally plunge into their central stars in less than a million years. In previous work, we have 
shown that a poorly ionized, less viscous region in a protoplanetary disk called a dead zone slows 
down the migration of fixed-mass planets. In this paper, we extend our numerical calculations to 
include dead zone evolution along with the disk, as well as planet formation via accretion of rocky 
and gaseous materials. Using our symplectic-integrator-gas dynamics code, we find that dead zones, 
even in evolving disks wherein planets grow by accretion as they migrate, still play a fundamental 
role in saving planetary systems. We demonstrate that Jovian planets form within 2.5 Myr for disks 
that are ten times more massive than a minimum mass solar nebula with an opacity reduction and 
without slowing down migration artificially. Our simulations indicate that protoplanetary disks with 
an initial mass comparable to the minimum mass solar nebula (MMSN) only produce Neptunian mass 
planets. We also find that planet migration does not help core accretion as much in the oligarchic 
planetesimal accretion scenario as it was expected in the runaway planetesimal accretion scenario. 
Therefore we expect that an opacity reduction (or some other mechanisms) is needed to solve the 
formation timescale problem even for migrating protoplanets, as long as we consider the oligarchic 
growth. We also point out a possible role of a dead zone in explaining long-lived, strongly accreting 
gas disks. 

Subject headings: accretion, accretion disks - turbulence - planetary systems: formation - planetary 
systems: protoplanetary disks - planets and satellites: general - Solar system: for- 
mation - stars: pre-main-sequence 



1. INTRODUCTION 

The properties of nearly 300 recently discovered extra- 
solar planetary systems reveal that Jovian mass planets 
are often found at a scaled orb ital radius of Mercur y 
around their central stars (e.g. lUdrv fc Santos! l2007t ) . 
Since none of the current theories of Jovian planet for- 
mation can explain in situ formation of gas giants at 
these small distances from their stars, it is generally 
agreed that such planets were formed in the outer re- 
gions of protoplanetary dis ks and migrated through them 
to their current positions (|Lin et al.lll996| ). Growth of 
giant planets takes place during this passage and both 
processes are terminated when most of the gas in the pro- 
toplanetary disk is either accreted onto the central star , 
or dissipated by photo e yaporation (e.g. Shu et all 119931 : 
iHollenbach et all 11994 I Alexander et alj l2006l ).~Obser- 
vations of infrared to submm emission and gas accre- 
tion rates onto the central stars reve al that disk life- 
times are typically 1 — 10 Myr (e.g. lHartmann et al.l 

1 This work was done while at McMaster University. 



1991: iMuzerolle et all [20001 : lAndrews fc Williams! [20051 : 
Sicil ia-Aguilar et al.ll2006 



The observed disk life-times raise two difficulties in 
planet formation theory. The first regards planet migra- 
tion: the migration time scales of planets arising from the 
tidal interaction between a planet and the gas disk are 
shorter than the disk lifetimes. Therefore, unless they are 
stopped by some robust mechanism, planets plunge into 
their central stars within about a million years. The sec- 
ond is specific to core-accretion model for Jovian planet 
formation: the formation time scales may be longer than, 
or comparable to, the disk life times. Therefore, unless 
the formation time scale is reduced somehow, the exis- 
tence of giant planets cannot be explained by the core 
accretion scenario. 

The core-accretion model posits that gas giants result 
from a two-stage process - the first being the formation 
of their rocky cores by the repetitive coagulation and ag- 
glomeration of the smaller bodies onto the larger ones (as 
in the terrestrial planet formation), and the second being 
the accretion of their massive gaseous envelope from the 
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surrounding disk. 

In the first stage, the formation of the rocky core pro- 
ceeds as bodies grow from a micron to a planetary size. 
One of the difficulties occurs when particles grow up 
to dm sizes. At this point, the collisional agglomera- 
tion m ay not be a preferred p ath to make the larger 
bodies (Eangkows ki et all [20 08). Even if the sticking is 
still efficient and the bodies can keep on growing, the 
gas drag becomes non-negligible for such objects. Since 
the gas disk rotates at the sub-Keplerian speed, slightly 
slower than these growing bodies, the bodies feel the 
head wind, lose the angular momentum, and eventu- 
ally migrate into the central star. Migration induced 
by gas drag becomes most efficient for meter-sized bod- 
ies, and becomes negligible again for the km-sized bodies 
since they are large enough no t to be affected by the 
gas drag (| Weidenschillind Tl977f ) . The gravitational in- 
stability in the planetes imal disk was suggested t o avoid 
this meter-size barrier (Goldrcich & Ward 1973). How- 
ever, such a mechanis m may be hindered by the Kelvin- 
Helmholtz instability (|Cuzzi et al.lll993t IWeidenschillind 
Il995f) unless the local solid-to-gas ratio is sufficiently high 
fe.g. lSekivalH998l : lYoudin fc ShiJl2f)02h . Recently, it was 
demonstrated that very rapid planetesimal growth up to 
Ceres-mass objects is possible via a streaming instability 
in such turbulent disks ( e.g. Uohansen fc You din 2007; 
lYoudin fe Joha nscn 2007). Therefore, there is good jus- 
tification to skip these early stages, and to assume that 
there are already km-sized planetesimals as well as plan- 
etary embryo(s) which are embedded in a gas disk. In 
this paper, we will follow this approach as the previous 
studies did. 

In the second stage, protoplanetary cores keep on ac- 
creting planetesimals while they start devel oping gaseous 
envelo pes. A cornerstone work done by iPollack et al.l 
(1996) (P96 hereafter) identified three phases in giant 
planet formation. The first phase is a rapid core building 
phase, in which a protoplanetary core of 0.6Me grows up 
to ~ 10Mb within half a million years. The second phase 
is a slow gas and planetesimal accretion phase which lasts 
until the crossover mass (for which the envelope mass is 
comparable to the core mass) is reached. The third phase 
is a rapid gas and planetesimal accretion phase, in which 
the planet quickly becomes a gas giant. P96 showed that 
planet formation spends most time in the second phase 
that could last for several million years. Since this is 
uncomfortably close to a typical disk life time, it was 
considered as one of the wea kest points of core accre- 
tion scenario (e.g. Boss] Il997t ). To shorten the second 
phase, protoplanets could either increase the gas accre- 
tion rate itself, or increase the planet esimal accret i on rat e 
and expand th e gas feeding zone. Ilkoma et all (|2000T) : 
iHubickvi etal] (|2005l) (hereafter INE00, and HBL05 re- 
spectively) took the former approach, and considered a 
reduced opacity in planetary envelopes. HBL05 used an 
updated version of the code by P96 with a more recent 
opacity table, and showed that smaller opacity, which 
mimics the dust settling and coagulation in the proto- 
planet's atmosphere, decreases the effective gas pressure 
of the protoplanetary envelope, and hence leads to the 
faster gas accretion. They showed that the in situ for- 
mation timescale of Jupiter could be as short as 1 Myr 
if the opacity is 2% of the interstell ar value and t h e core 
mass is 10M E . On the other hand. lAlibert et alJ (|2005l ) 



(hereafter A05) took the latter approach, and considered 
planet formation during migration. In P96, planetesimal 
accretion slows when the planetesimals get depleted in 
the feeding zone, because a planet has to accrete gas and 
expand its feeding zone to further accrete planetesimals. 
Since the gas accretion in Phase 2 tends to be slower than 
the planetesimal accretion in Phase 1, it takes longer to 
accrete a similar amount of gas envelope to the core, 
and achieve a crossover mass. A05 overcame this prob- 
lem by including disk evolution and planet migration in 
the model, so that the planetesimal feeding zone is con- 
stantly replenished. The gas accretion time shortens as 
the planetary mass increases, because it proceeds on the 
Kelvin-Hclmholtz timescale, which is a steep function of 
planetary mass. They showed that the Jupiter formation 
timescale could be of the order of 1 Myr if t he planet mi- 
grates. However, as noted by IWardl (|1997l ). solving the 
accretion problem by migration in a standard disk model 
is a "double-edged sword" , since it comes at the price of 
possibly losing the planets to their central star. Thus, 
A05 assumed that some unspecified mechanism would 
slow down the fast type I (pre-gap opening) migration 
by factors of 10 — 100 tim es compared to the migration 
estimated in a 3D disk bv lTanakaet~aT1 (|2002t ). 

Here, we propose a possible mechanism which can slow 
planet migration, and investigate planet formation in 
that context. In a standard disk with a smooth surface 
mass density and a standard viscosity expected from the 
MRI turbulence, both type I and II migration time scales 
are shorter than the disks' life times. The most popu- 
lar source of such a viscosity is the magn eto-rotational 
instability (MRI. iBalbus fc Hawlevl 119911 ). However, it 
has been demonstrated by several groups that proto- 
planetary disks are not MRI active everywhere, but har- 
bor e xtended regions known as the dead zones (|Gammiei 
Il996| ) where there is virtually no turbulence 2 . Many au- 
thors have studied the physical extent of the dead zones 
by using different disk models and taking account of a 
variety of the ionizat i on/recombination processes (e.g. 
iGlassgold et alJ [199ft ISano et all 12000] lEYomang et al l 
| 2002t ISemenov et al.ll2"00l iMatsumura fc Pudritzi l200l 
I2OO6I hereafter MP03, and MP06 respectively). These 
models suggest roughly the same size of the dead zones 
(from < 1 AU to 10 — 20 AU), indicating that a dead zone 
is a robust feature of a protoplanetary disk. Since this is 
a critical region of a disk both for planet formation and 
migration, it is of vital importance to investigate these 
problems in the context of protoplanetary disks with 
dead zones. Recent numerical simulations of such a lay- 
ered disk have shown that the expected viscosity param- 
eter a is about 10 -2 in the active zone and 10~ 4 — 10~ 5 in 

2 This arises in principal because operation of the MRI requires 
good coupling between the gaseous disk and the magnetic field. 
However, the ionization fraction in the inner regions of dense pro- 
toplanetary disks is typically so low that Ohmic diffusion prevents 
the growth of the MRI in such regions. The size of the dead zone 
can be computed once the source of disk ionization (e.g. stellar 
X-rays, cosmic rays) is known. Assuming that the MRI is the 
dominant source of the disk's "turbulent" viscosity, the dead zone 
is nearly inviscid. The existenc e of a de a d zon e in a protoplane- 
tary disk was first proposed by Gammic (1996) who showed that 
there is a magnetically dead zone in the inner disk where cosmic 
rays cannot penetrate (assuming the cosmic ray stopping density 
is 98 g cm" 2 ), and that such a region is sandwiched by the mag- 
netically active surface layers where gas accretes toward the star 
efficiently. 
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the dead zone (jFleming fc Stondl2003f) . Therefore, if the 
MRI is the dominant source of the disk viscosity, a sig- 
nificant decrease in viscosity is expected within the dead 
zone, which most likely affects the rate of planet migra- 
tion. Moreover, recent observations have revealed that 
disk's viscosity parameters may take a larger range than 
previously expected: a = 10~ 6 — 10 _1 (|Hueso fc Guillotl 
l2005t iLuhman et al.ll2007t ). indicating that at least some 
observed disks may possess a low viscosity region like a 
dead zone. 

In an earlier paper, we have demonstrated that dead 
zones can significantly slow the rate of migraion of plan- 
ets, and save planetary systems from plunging into the 
central stars (|Matsumura et al.ll2007l hereafter MPT07). 
In particular, we found (1) type II migration (post-gap- 
opening migration) is slowed in the dead zone due to the 
low viscosity there, (2) even low-mass planets (< lOMg), 
which are usually type I migrators (i.e. non gap-openers), 
may open a gap in the dead zone if the thermal condition 
is satisfied, and thus migrate slower there, and (3) type I 
migrators moving toward the dead zone can be stopped 
at the outer edge of the dead zone due to the jump in 
mass density there which is a result of the slower advec- 
tion speed inside the dead zone with respect to outside 
it. 

In this paper, we present a rather complete treatment 
of the formation and migration of Jovian planets within 
their evolving protoplanetary disks. We compute com- 
prehensive time-dependent models for the growth of Jo- 
vian planets in the core accretion picture in viscously 
evolving gaseous disks. Our calculations include updates 
on both stages of accretion (planestimals and gas), as 
well as the ability to follow planetary migration down to 
0.1 AU. We approach this problem by means of time- 
dependent simulations that track both planetary accre- 
tion and migration through disks with evolving dead 
zones. There are two major effects of disk evolution in 
our models: (i) the gradual accretion of mass from the 
disk onto the star which limits the reservoir that is avail- 
able to the growing Jovian planet; and (ii) the shrinkage 
of the size of the dead zone as a consequence of mass 
accretion onto the star (mainly) through the well cou- 
pled surface layers of the disk. This process reduces the 
size of the dead zone substantially with time — perhaps 
leaving it only a few AU in extent after a million years. 

We first introduce our disk models and numerical 
methods in §2. We also highlight a possible role of dead 
zones in disk accretion onto the central stars. Then we 
study planet formation in an inviscid disk, and compare 
our results with P96, and HBL05 (§3). We also check 
the effect of the opacity, and compare the results with 
HBL05, which improved the work of P96 and further in- 
vestigated the opacity effects. We then go on to compute 
planet formation and migration in an evolving disk, and 
compare the results with A05 in §4. The culmination 
of our work is presented in §5 where we generalize our 
results and present planet formation in an evolving disk 
with a dead zone. Finally, we summarize our work in §6. 

2. NUMERICAL METHODS AND DISK MODELS 

In this section, we introduce the numerical code that 
we use to simulate planet formation in a realistic proto- 
planetary disk. We focus on the newly added features 
and refer the reader to our previous paper (MPT07) for 



details of the basic code. We also summarize the initial 
conditions chosen for the runs in the following sections. 

We perform numerical simulations of planet forma- 
tion and migration by using a hybrid numerical code, 
which combines an N-body integra tor with a simple 
disk evolution code (Thommes 2005). The N-body part 
is based on the sy mplectic integrator called SyMBA 
(jDuncan et alj I1998D. which has im proved the N-body 
map of IWisdom fc Holmanl (fl99lh with an adaptive 
timestep to handle the close encounters among massive 
bodies. The gaseous disk part of the code evolves a 
disk viscously as well as through angular momentum ex- 
change with the embedded planets according to a gen- 
eral Navier-Sto kes equation. By f ollowi ng the standard 
prescription bv lLin fe Papaloizoul (|1986t ). the gas disk is 
divided into radial bins, which represent annuli of disks 
with azimuthally and vertically averaged properties like 
surface mass density, temperature, and viscosity. Viscous 
evolution of the disk is calculated by specifying the stan- 
dard alpha viscosity (jShakura fc Sunvaevlll973D for each 
bin, while the effect of disk-planet interac tions i s adde d 
in the form of the t orque density as in IWardl (fl997h : 
IMenou fc Goodman! (|200l . 

The code is also modified to simulate planet forma- 
tion, which consists of planetesimal and gas accretion. 
For planetesimal accretion, as it was mentioned in §2, the 
transition between runaway, and oligarchic phases occurs 
once the protoplanet becomes locally large enough to per- 
turb the surrounding smaller planetcsimals and increase 
their random velocities, which decreases the collisional 
crosss ection, and therefore leads to the longer a ccretion 
time (|Ida fc Makind fl99l IKokubo fcldal fl99l . This 
criti cal mass is about a few times 10~ 3 Moon masses or 
less |Th ommcs "etaLl [2003 s ). Since all of our simulations 
start with 0.6 Me, we can safely assume that the oli- 
garchic growth is the dominant planetesimal accretion 
phase in our case. 

The plan etesimal accretio n part of the code was de- 
veloped by iThommes et al.l (|2003h to handle the oli- 
garchic growth, which is the slower planetesimal accre- 
tion phase following the rapid planetesim al accretion 
(|Ida fc Makind 11993 IKokubo fc Idal Il998h. Her e, we 
follow the previous s tudies (e.g. IKokubo fc Idal 119961 : 
IThommes et al.1 12003T ) . and assume (1) dynamical fric- 
tion by smaller planetesimals on larger ones is effective, 
so that larger planetesimals have smaller velocities than 
smaller ones, and (2) gravitational focusing is effective, 
so that the relative velocities between two colliding bod- 
ies are smaller than the escape velocities from the larger 
ones, and therefore collisions effectively lead to coagula- 
tion. Accretion rate of planetesimals with mass m onto 
a protoplanet with mass M p is written as 



dM p 
dt 



solid-Mp 



4/3 



,1/2 



(1) 



where C = 67r 2 / 3 [3/(4p M )] 1/3 (G/M*) 1/2 , pM i s the bulk 
density of a protoplanet, I^ so iid is the surface mass den- 
sity of a planetesimal disk, and e m is the eccentricity 
of planetesimals. We assume that the protoplanet ac- 
cretes planetesimals within its feeding zone of 10i?Hiii 
(Rmu = (M p /(3M 4 )) 1 / 3 fl is the Hill radius of the pro- 
toplanet), whi ch is a typical orbit al separation between 
protoplanets (jKokubo fc Idalll998l ). 
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The gas accretion part of the code estimates the plan- 
etary mass increase by calculating the gas accretion 
timescale and radius (see §2.1 for details). Instead of 
solving the planetary structure equations directly to es- 
timate the gas accretion rate as in P96, HBL05, and 
A05, we calculate the gas accretion timescale based on 
these previous studies, and let the protoplanetary core 
accrete gas inside the accretion radius on the estimated 
timescale. 

We also add a new capability to our disk model to 
include the evolution of a dead zone due to the faster 
accretion through the well coupled surface layers onto 
the central star. This is further explained in §2.2. 

2.1. Gas accretion prescription 

Since the planetesimal accretion prescription is dis- 
cussed in detail elsewhere (jThommes et al.ll2003l ). we fo- 
cus on the gas accretion prescription added to the code 
in this subsection. Gas accretion onto a protoplanet 
has been studied either by solving the planetary struc- 
ture equations (e.g. P96, HBL05, INEOO, and A05), 
or b y performing the hydrodynamic (HP) simulations 
(e.g. iBrvden et alJ Il999l iLubow et all 119991: iKlev et all 
20011: ID'Angelo et alJl2002t iT anigaw afe Watanabdl2002t 
D'Angelo et alj|2003t iBate et alN2003fh The main focus 
of the former studies is on the core-building phase to the 
rapid gas accretion phase, while that of the latter ones is 
on the final stage of gas accretion where a massive planet 
accretes from its subdisk. Since the former studies ex- 
cept A05 treat the gas disks rather simply by assuming 
that the protoplanetary feeding zone is always replen- 
ished, they have a relatively crude estimate of the avail- 
able amount of gas in the feeding zone. Therefore, they 
tend to overestimate the gas accretion for a planet which 
is massive enough to open a gap. The gas accretion in 
the latter case is regulated by the subdisk formed around 
the planet. Since these studies do not include the con- 
traction of the protoplanetary atmosphere, they tend to 
overestimate how quickly a protoplanetary core can ac- 
tually accrete gas. Here, we try to combine the strengths 
of these two methods, and use three different gas accre- 
tion rates depending on presence/ absence of planetesimal 
accretion, and circumplanetary disk accretion. 

The presence, or absence of planetesimal accretion af- 
fects the efficiency of gas accretion. INEOO suggested 
that, when the planetesimals are accreting concurrently 
with the gas, the thermal energy released in the bom- 
bardment of planetesimals provides extra pressure sup- 
port for the gaseous envelope of the protoplanet, and 
therefore the gas accretion slows down. They found that 
such an accretion rate could differ from the accretion rate 
without planetesimal accretion by several times, which 
are plotted in Fig. Q] as two parallel lines. They also 
studied the effects of core mass and dust opacity on 
gas accretion, and showed that there is an optimal core 
mass which results in efficient gas accretion. If plan- 
etesimal accretion is cut off long after the protoplanet 
starts accreting gas actively, the gas accretion time be- 
comes longer. This is because the radiative loss is mostly 
compensated by planetesimal accretion rather than gas 
accretion (i.e. gravitational energy release due to the en- 
velope contraction). On the other hand, if planetesimal 
accretion is cut off while the core is still small, the gas ac- 
cretion time is also prolonged, because the smaller mass 



leads to the smaller radiative loss and hence the slower 
envelope contraction. 

The final stage of gas accretion is likely to be controlled 
by th e subdisk around a planet dTanigawa fc Watanabd 
120021 : IBate et al.1 I2003L ID'Angelo et all 120031 ). Gas 
accretion through subdisks becomes less efficient for 
more massive planets b ecause the planet s open a wider 
and deeper gap (e.g. IBate et all [2003). The effect 
is seen in the para bola-like curve in Fig. [T] by 
ID'Angelo et al.l (120031). This , however, is not clear in 
(jTanigawa fc Watanabd l2002t hereafter TW02)'s case 
(the lower dotted line), because they don't take account 
of a gap-opening effect. In other words, their gas accre- 
tion rate is the "raw" accretion rate of the planet. Be- 
low the critical mass for subdisk accretion, gas accretion 
is expected to be controlled by the contraction of pro- 
toplanetary atmosphere. Therefore, we apply the KH 
timescale for small protoplanets, and subdisk accretion 
timescale for larger protoplanets. 

Following these studies, we define three stages of the 
gas accretion; (1) gas accretion while planetesimal ac- 
cretion is on-going, (2) gas accretion after planetesimal 
accretion ceases, and (3) gas accretion through the sub- 
disk around the planet. Fig. [T] presents a compilation 
of these gas accretion timescales obtained from differ- 
ent studies. Two parallel descending lines are gas accre- 
tion time scales while planetesimal accretion is on-going 
(upper line) and after it ceases (lower line) respectively 
(INEOO, see Eq. (1) and (2)). For comparison, accre- 
tion time scales obtained from P96 is shown in the upper 
dotted line. These descending trends predict shorter ac- 
cretion timescales for larger planets because the models 
don't take account of disk evolution effects nor planetary 
subdisks. On the other hand, the parabola- like re l ation i s 
based on the HD simulation by ID'Angelo et alj (2003) , 
which predicts much longer gas accretion timescale for 
larger planets due to gap opening and accretion via cir- 
cumplanetary disks. A similar study done by TW02 
without a gap-opening effect is shown in the lower dotted 
line. Possible gas accretion timescales throughout planet 
formation are shown in solid curves. 

(1) When there is planetesimal accretion as well as gas 
accretion, we follow the timescale by INEOO (the upper 
line of the solid parallel lines): 

-2.5 



n = 6. x 10 8 



M r , 



(2) 



M*J \1 cm 2 g" 
where k is the dust opacity in the planetary envelope. 

(2) When the planetesimal accretion ceases, but the 
gas accretion is still on-going, we adopt the timescale by 
INEOO (the lower line of the solid parallel lines): 

— ° 8 (sHot)- < 3 » 

Both of the above equations show shorter accretion 
timescales as the planetary mass grows. 

(3) When a protoplanet becomes large enough to have 
a s ubdisk, the gas accret ion timescale can be described 
bv lD'AngeloetaH (l2003h : 



T3 - 



Mp 



log 



M v 



M E yr- 



18.47 + 9.25 log ( ^ 



1.2661og g 
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The transition of the first to the second phase occurs 
when the core accretion rate drops to zero. The transi- 
tion of the second to the third phase occurs when ti ~ T3. 
For a standard opacity k ~ 1 cm 2 g _1 , the crossover 
mass is about 70 Me as it can be seen in Fig. [TJ 

A protoplanet accretes gas within its accretion radius 
on each of these time scales. As in P96, the accretion 
radius for Eq. [2] and [3] is set to either Hill (Rmu = 
(Af p /(3A/»)) 1/3 a), or Bondi (i? B ondi = GA/»/c 2 , where 
c s is sound speed) radius, whichever is smaller. On the 
other hand, the accretion radius for Eq. [5] is chosen to 
be 2-Rniii, which is m otivated by numerical studies like 
iD'Angelo et "all (|200l . 

2.2. Accretion histories in disks with dead zones 

Although our disk model is one-dimensional, it is still 
possible to include the effects of azimuthal and vertical 
structure through disk parameters in an averaged way. 
The difference in disk evolution due to a vertical struc- 
ture becomes particularly important for a disk with a 
dead zone. Generally, a dead zone is sandwiched between 
the upper and lowe r turbulent sur face layers where the 
disk is well- ionized (|Gammidll996l) and hence the MRI 
is active. Therefore, the mass accretion is expected to be 
more efficient in active layers than in the dead zone. In 
MPT07, we did not include this effect, and hence found 
very little decrease in disk mass even after 10 7 years. 

Here, we take account of the effects of disk's verti- 
cal structure in the following simple way. The mass 
accretion through the active layers and the dead zone 
toward the central star can be written as the summa- 
tion of accretion rates in the active layers and the dead 
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Fig. 1. — Comparison of gas accretion timescales as a function 
of planetary mass. The descending solid and the upper dotted 
line are the gas accretion timescales for the envelope contraction 
phase, while a parabola and the lower dotted line are those for the 
subdisk accretion phase. The upper dotted line is the gas accretion 
rate estimated from the simulation of P96. Blue line (the upper 
of the parallel lines) is gas accretion while planetesimal accretion 
is on-going, while black one (the lower of the parallel line) is that 
after planetesimal accretion ceases. Both of them are based on 
the study by INE00. Black parabola an d the lower dott e d line 
are the subdisk accretion timescales from D'Angclo ct al. (2003) 
and T W02, respectively. The subdisk accretion by D'Angclo ct al. 
ll2003l) appears to be slower than TW02 since the gap-opening effect 
of a planet is included. 
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Fig. 2.— The surface mass density at t = 0, 10 4 , 10 s , 10 6 , and 
10 7 years for disk with dead zone (DZ) but no planet (RunA). The 
dead zone shrinks as the disk evolves. 

zone (e.g. ILvnden-Bell & Pringld H97l IGammid [1996; 
iFromang et al.ll2002f ): 

= 67rr 1/2 ^-(2S crit ^active?" 1/2 +(£-2i: cr it) ^doad? ,1/2 ) , 
at or 

where £ C Ht is the critical surface mass density below 
which the disk is well-ionized, while £ is the total gas 
surface mass density. Specifically, we assume that the 
local disk is "dead" if £ > £ cr ;t = 21, and 80 g cm~ 2 
for the surface mass density at 1 AU of So = 10 3 and 
10 4 g cm -2 , respectively. The factor of 2 in front of £ cr it 
implies the upper and lower layers. Also, v = ac s h is the 
viscosity (a is the viscosity parameter, c s is the sound 
speed, and h is the disk pressure scale height), and the 
subscripts dead and active imply parameters in the dead 
zone and active layers respectively. 

From the above equation, we can define the "vertically 
averaged" disk viscosity parameter that is weighted by 
the surface mass density: 

a = £ ( 22 "' critaactivc + ( S _ 2£ crit ) adcad) . (7) 

In our simulations, the protoplanetary disks evolve vis- 
cously according to this viscous a. 

Now we will show one example which describes a typ- 
ical evolution of a disk with a dead zone. Using the 
averaged alpha value in Eq. [3 and assuming that 
a a ctivo = 10~ 2 and «dcad = 10 -5 , we evolve a MMSN- 
type disk with £ = 10 3 (a/AU)~ 3/2 gem" 2 as in Fig. H 
and [3] (RunA, see Table [1] for initial conditions). Fig. 
[2] shows the surface mass density evolution of a disk with 
a dead zone, including the effect of the surface layer ac- 
cretion. There is no planet in this case. An important 
feature of disks with dead zones is immediately apparent 
in this figure. The high advection speed in the outer disk, 
compared to the much lower speed in the region contain- 
ing the dead zone, results in a pile up of gas. This man- 
ifests itself in a steep density gradient, which we have 
already shown can play an important role in reflecting 
the inwardly migrating planets in the outer disk regions 
(MPT07). As material accretes along the surface of the 
dead zone, we see that this dense front moves inward. 
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Fig. 3. — Evolution of disk with DZ but no planet (RunA) Top: 
Evolution of the disk column density. The disk viscosity parameter 
is set to a = 10 -5 inside the dead zone, and 10 — 2 outside it. The 
dead zone shrinks rapidly due to mass accretion through the surface 
layers. Bottom: Corresponding plots of disk accretion rate (solid 
line) and disk mass (dashed line). Once the dead zone is gone, the 
disk disappears quickly. 



The dead zone shrinks very rapidly, and only within 2 
AU is dead by 2 Myr. We will see how the evolution of 
a dead zone affects planet formation in §5. 

The corresponding figure for the disk evolution is 
shown in the top panel of Fig. [3J Note that the dead 
zone edge is where the surface mass density contours are 
the densest. The accretion onto the central star speeds 
up as the dead zone disappears. The evolution of mass 
accretion rate (solid line) and disk mass (dashed line) 
for the same run are shown in the bottom panel of Fig. 
[3J With a dead zone, the mass accretion rate is nearly 
constant for a few Myr, and more than 10% of the initial 
disk mass accretes onto the central star once the dead 



zone is gone at around 10 Myr. This final dispersal of 
a disk proceeds at the viscous time scale of a turbulent 
disk if there is no other mechanism, and takes about a 
few Myr. 

This characteristic evolution of the disk accretion rate 
in disks with dead zones may have an interesting im- 
plication for the observations. The current observa- 
tions indicate that mass accretion rate is a gradually 
decre asing function in time (e.g. iSicilia-Aguilar et "all 
2006), although the range in mass accretion rate 10 -10 — 
10~ 7 M Q yr~ 1 is roughly constant between 1 — 10 Myr. 
In a standard disk with no dead zone, a gradual decrease 
of disk mass is expected (see Fig. [B]for example), and 
therefore we may need very massive initial disks to ex- 
plain the long-lived, strong accretors. 

Alternatively, if the dead zone is playing a role in disk 
dispersal, the mass accretion rate is likely to stay nearly 
constant throughout the evolution until the dead zone 
disappears. In this case, mass accretion rates correlate 
with disk masses as well as the difference between effec- 
tive viscosities inside and outside the dead zones. Once 
the dead zone disappears, the rest of the disk accretes 
onto the central star at a higher viscous accretion rate. 
Currently, it is difficult to say whether gas accretion rates 
tend to decrease gradually, or stay rather constant for a 
significant part of the disks' lifetime. Future observations 
may constrain masses and accretion rates of protoplanc- 
tary disks more precisely to indicate a possible preferred 
path. 

2.3. Disk models 

We study planet formation in a protoplanetary disk by 
assuming that a single planetary embryo with 0.6 Me is 
embedded in gas and planetesimal disks. This is the 
same initial core mass used in P96, HBL05, and A05. In 
this paper, we use three slightly different gas disk models 
to compare our results directly with the previous studies, 
but keep the solid surface mass density the same S so iici = 
270(r/AU)~ 2 for all the models. Instead of choosing a 
specific planetesimal size, we use different planetesimal 
sizes ranging from 100 m to 100 km, and compare their 
results with one another. 

In §3, we check our gas accretion prescription against 
P96 and HBL05, and simulate the formation of a planet 
on a fixed orbit in an inviscid disk. To make the 
comparison easier, we adopt the same initial condi- 
tions as P96 and HBL05 — the viscosity parameter is 
a = 0, the gas surface mass density profile is E gas = 
700(r/(5.2 AU))- 2 = l8950(r/AU)- 2 , the disk temper- 
ature is T = 150(r/5.2AU)-^ 2 = 342(r / AU)- 1 ' 2 , and 
the disk extends from 0.25 to 50 AU. We use both a stan- 
dard k = 1 cm 2 g -1 and a reduced k = 0.03 cm 2 g _1 
opacity for our runs. 

In §4, we include the effect of disk evolution and planet 
migration, and compare our results with A05. Again, we 
choose the same initial conditions as A05 to make the 
comparison easier, which is the same as §3 except a = 2x 
10~ 3 , and £ flas = 525(r/(5.2 AU))- 2 = Ul96(r/AU)- 2 . 

Finally in §5, we study the effect of a dead zone on 
planet migration and formation. We use the same disk 
model as MPT07 — E gas = Y. (r/AU)- 3 / 2 with S = 
10 3 or 10 4 g cm" 2 , an d the disk temperature model by 
iRobberto et alj (|2002f ) , which takes account of the effect 
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of a cluster environment. In this model, the disk extends 
from 0.02 to 100 AU. Throughout this section, we use 
the viscosity parameter a = 10~ 5 inside the dead zones, 
and a = 10~ 2 for the disk beyond it (active zones). Here, 
the initial value of the outer dead zone radius is ~ 8.2, 
and 15.8 AU for Eo = 10 3 , and 10 4 g cm -2 , respectively. 
These values of viscosity para meter agree well with th e 
numerical simulations done bv lFleming fc Stoni (|2003h . 

The initial conditions of all runs are summarized in 
Table [U 

3. PLANET FORMATION IN NON-EVOLVING 
PROTOPLANETARY DISKS 

In this section, we check our gas accretion prescription 
against previous studies. We simulate in situ planet for- 
mation in a non-evolving protoplanetary disk by adopt- 
ing the same initial conditions as P96 and HBL05, and 
compare our results with theirs. The disk models in this 
section assume an inviscid disk (i.e. disk viscosity pa- 
rameter is a = throughout the disk) which extends 
from 0.25 to 50 AU (see §2.3 for details). 

3.1. Effect of the planetesimal size 

In this subsection, we focus on the effect of the plan- 
etesimal size on planet formation timescale. Generally 
speaking, smaller planetesimals lead to faster accretion 
because they are subject to the stronger damping of ran- 
dom velocities, and hence have larger cross-sections for 
accretion by a protoplanet. 

P96 investigated the giant planet formation in a non- 
evolving disk, and found that a Jupiter-like planet can 
form in situ within about 8 Myr, while a similar study 
with an improved equation of state and opacity tables by 
HBL05 predicts 6 Myr. We use the same initial condi- 
tions as P96 and HBL05, and run four different simula- 
tions by changing only the planetesimal size from 100 m 
to 100 km (RunBl-4). RunB4 corresponds to the fiducial 
model in P96, and lOtfoo in HBL05. 

Fig. d] shows the results of our simulations of gas 
and planetesimal accretion by a protoplanetary core with 
0.6Me at a fixed orbit (5.2 AU). The top panel is the time 
evolution of core, envelope, and total masses of a proto- 
planet for a planetesimal size of 100 m, 1 km, 10 km, 
and 100 km from left to right. The bottom panel shows 
the time evolution of core and envelope mass accretion 
rates for the corresponding runs. Once the core accre- 
tion ceases, the transition from Eq. [2] to [3] occurs, and 
the envelope accretion rate increases, which is marked by 
a jump. Planets approach their final masses as the gas 
accretion slows down due to the subdisk accretion and a 
gap-opening. 

As the previous studies have shown, we find that the 
core building takes longer for a swarm of larger planetes- 
imals, and so does the time until the start of the rapid 
gas accretion. As it was mentioned earlier, this is be- 
cause smaller planetesimals have larger accretion cross- 
sections by a protoplanet. For 100 m-size planetesimals, 
a Jupiter-like planet forms within about 2 Myr, while for 
10 km-size planetesimals, it takes about 7 Myr, which 
is roughly the same as the time estimated for the fidu- 
cial model of P96 (~ 8 Myr) or the corresponding model 
(lOifoc) in HBL05 (~ 6 Myr). 

Now we compare our results with the fiducial model 
of P96 and lOiioo in HBL05. P96(HBL05) obtained the 
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Fig. 4. — Evolution of planetary masses and accretion rates for 
planet at fixed disk orbital radius of 5.2 AU (RunBl-4). Top: Evo- 
lution of planetary mass which is initially O.GMe- The planet has 
a fixed orbit at 5.2 AU from the central star. The black curves 
show the total masses of planets, the orange curves show the core 
masses, and the blue curves show the envelope masses. Different 
sets of curves represent planetary growth with different planetesi- 
mal sizes (100 m, 1, 10, and 100 km.) The accretion time becomes 
shorter for smaller planetesimals. The models with planetesimal 
size about 10 km produce a similar result to P96, and 10-ffoo in 
HBL05. Bottom: Corresponding mass accretion rates. Orange 
curves are core accretion rates, and blue curves are gas accretion 
rates from Eq. |2][3] and [5] 

estimated formation timescale of 8(6) Myr for 100 km- 
size planetesimals. In our run, the corresponding case 
(RunB4) did not achieve a Jupiter mass within the sim- 
ulation time (10 Myr). Moreover, although the mass 
evolution profiles look qualitatively similar to the ones in 
P96 or HBL05, the details appear to be different. Phase 
1 and 2 in their cases are merged in our case. The differ- 
ence is apparent by comparing the bottom panel of Fig. 
[4] with Fig. lb in P96. In their case, rapid planetesi- 
mal accretion (Phase 1) ceases when a protoplanet de- 
pletes its planetesimal feeding zone, and this is followed 
by slow gas accretion (Phase 2), which lasts until the 
planet obtains the crossover mass for rapid gas accretion 
(Phase 3). Since they assume rapid planetesimal accre- 



TABLE 1 

Initial conditions for all the runs. From left to right columns, there are run name, initial gas surface mass density 
(S gas = E(j(a/AU)~ p [g cm" 2 ]), DISK temperature (T = To(a/AU)~ q [K]), planetesimal size, viscous alpha values in active and 

DEAD ZONES, OPACITY, AND TORQUE SOFTENING PARAMETER. 



Run 


(S [gem- 2 ], p) 


(To [K], ,) 


planeteismal size 


(^active, C*d cad) 


re [cm 2 g 1 ] 


B 


A 


(10 3 , 1.5) 


MPT07 


N/A 


(10~ 2 , 10~ b ) 


N/A 


N/A 


Bl 


(1.9 x 10 4 , 2) 


(342, 0.5) 


100 m 


(0, N/A) 


1 


N/A 


B2 


(1.9 x 10 4 , 2) 


(342, 0.5) 


1 km 


(0, N/A) 


1 


N/A 


B3 


(1.9 x 10 4 , 2) 


(342, 0.5) 


10 km 


(0, N/A) 


1 


N/A 


B4 


(1.9 x 10 4 , 2) 


(342, 0.5) 


100 km 


(0, N/A) 


1 


N/A 


CI 


(1.9 x 10 4 , 2) 


(342, 0.5) 


100 m 


(0, N/A) 


0.03 


N/A 


C2 


(1.9 x 10 4 , 2) 


(342, 0.5) 


1 km 


(0, N/A) 


0.03 


N/A 


C3 


(1.9 x 10 4 , 2) 


(342, 0.5) 


10 km 


(0, N/A) 


0.03 


N/A 


C4 


(1.9 x 10 4 , 2) 


(342, 0.5) 


100 km 


(0, N/A) 


0.03 


N/A 


Dl 


(1.4 x 10 4 , 2) 


(342, 0.5) 


10 km 


(2 x 10- 3 , N/A) 


1 


0.6 


D2 


(1.4 x 10 4 , 2) 


(342, 0.5) 


100 m 


(2 x 10~ 3 , N/A) 


1 


0.6 


D3 


(1.4 x 10 4 , 2) 


(342, 0.5) 


100 m 


(2 x 10~ 3 , N/A) 


1 


0.9 


El 


(10 3 , 1.5) 


MPT07 


100 m 


(10~ 2 , 10~ 5 ) 


1 


0.6 


E2 


(10 3 , 1.5) 


MPT07 


100 m 


(10- 2 , IO" 5 ) 


0.03 


0.6 


Fl 


(10 4 , 1.5) 


MPT07 


100 m 


(10~ 2 , 10~ b ) 


1 


0.6 


F2 


(10 4 , 1.5) 


MPT07 


100 m 


(10~*, IO -5 ) 


0.03 


0.6 


Gl, & G2 


(10 3 , 1.5) 


MPT07 


100 m 


(10~ 2 , IO -5 ) 


0.03 


0.6 


HI 


(10 4 , 1.5) 


MPT07 


100 m 


(io- 2 , io- 5 ) 


0.03 


0.6 


H2 


(10 4 , 1.5) 


MPT07 


100 m 


(10~ 2 , IO" 5 ) 


1 


0.6 



tion, the planetesimal accretion rate is very high during 
Phase 1 (dM/dt ~ 1CT 4 - 1(T 2 M E /yr), and drops to 
~ 1CP 6 MEyr" 1 during Phase 2 due to planetesimal de- 
pletion in the feeding zone (see Fig. lb of P96.) Their 
protoplanetary core grows from 0.6Me to ~ lOAfg dur- 
ing Phase 1 (within 0.5 Myr). 

However, as we mention in §2, protoplanets more mas- 
sive than a few times 10~ 5 Me are expected to grow oli- 
garchically. In our simulations, since we assume slower, 
oligarchic growth from the start, the planetesimal accre- 
tion rate is low all the time ~ 1CP 5 — 10~ 6 Msyr -1 and 
drops to zero once the planetesimals are depleted. As 
a result, the core mass does not reach ~ lOAfg for a 
few Myr in RunB3. The difference in core growth pre- 
scription may not seem to have a significant effect on 
a non-migrating planet in the inviscid disk, but becomes 
important for a migrating planet in an evolving disk. We 
discuss this further in §4 and 5. 

3.2. Effect of the opacity 

A reduced opacity in the gaseous envelope of a pro- 
toplanet can lead to the faster gas accretion (P96 and 
INE00). HBL05 extended the work of P96, and stud- 
ied the effects of disk's opacity as well as a planetary 
core mass. By assuming 2% of the opacity of interstellar 
grains, they found that the gas accretion time becomes 
significantly shorter (1 — 4.5 Myr). This effect is reflected 
in the gas accretion timescale we adopt (see Eq. [2] and 

ED. 

We keep the definition of our solid surface mass den- 
sity, and change the size of planetesimals from 100 m 
to 100 km as in the last subsection (RunCl-4). RunC4 
corresponds to the model of lOLoo in HBL05. Here, we 
adopt a reduced, constant opacity of k = 0.03 cm 2 g -1 , 
which is an average opacity value obtained for the model 
10-Loo in HBL05 over a relevant range of temperature of 
protoplanetary disks (roughly 100 — 1000 K). 

The results of our simulations are shown in Fig. 
Comparing Fig. [?]with[ni we confirm that the lower opac- 
ity results in the faster gas accretion. HBL05 showed that 



the formation time could be as short as 2 Myr in a disk 
with 100 km planetesimals with a reduction in opacity 
(see their Fig. 1). In our simulations, a similar time scale 
to their model is obtained for a 10 km-size planetesimal 
disk (RunC3), where formation takes about 3 Myr. This 
is more than twice faster than our corresponding run in 
§3.1 with the standard opacity (RunB3). The formation 
time scales for our 10 km planetesimal disk models are 
about 1 Myr longer than those of corresponding models 
in HBL05 with a 100 km-size planetesimal disk (lOLoo 
and 10-ffoo). Thus, our gas accretion prescription agrees 
fairly well with previous studies like P96 and HBL05. 

4. PLANET FORMATION IN AN EVOLVING DISK 

We now add disk evolution and planet migration to the 
formation model described in the previous section, and 
study their effects on planet formation to compare the 
results with A05. In this section of the paper, the disk 
viscosity parameter is assumed to be constant (a = 2 x 
10~ 3 ) throughout the disk. Therefore, different from the 
previous section, a protoplanetary disk accretes toward 
the central star. Note that this model does not include 
a dead zone. The other initial conditions of the disk 
models are explained in §2.3, and summarized in Table 

ED 

A05 adopted th e migration speed obtained by 
iTanaka et al.l (|2002t ). and reduced the migration rate ar- 
tificially by a factor of 0.1 — 0.01. To incorporate the idea, 
but to keep the generality, we use the t orque expression 
presented bv lMenou &: Goodman! (|2004[ ). and adjust the 
softening parameter B, which mimics the torque reduc- 
tion effect with a vertical height (also see Appendix B in 
MPT07). 

First, we show our fiducial migration case with B = 
0.6, which ro ughly corresponds to the migration speed 
estimated by ITanaka et al.l (|2002t ). Fig. [5] shows the 
growth of a planet with an initial mass of 0.6Me in a 10 
km-size planetesimal disk (RunDl). The planetesimal 
accretion is very slow as it can be seen in the bottom 
panel, and therefore the protoplanet does not migrate 
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Fig. 5. — Effect of reduced opaciy on evolution of planetary 
masses and accretion rates for planet at fixed orbital radius of 5.2 
AU (RunCl-4). Top: The same as Fig. [4] but with a disk opacity 
of k = 0.03 cm 2 g — 1 . The model with planetesimal size about 10 
km results in a similar outcome to the model lOLoo in HBL05. 
Bottom: Corresponding mass accretion rates. 

much. The final mass is about 1 Me, only barely in- 
creased from 0.6 Me- Also plotted in the bottom panel 
of Fig. [6] are the disk mass evolution (dashed line) and 
disk accretion rate onto the central star (solid thin line) . 
It is interesting to compare this figure with the disk evo- 
lution with a dead zone (see the bottom panel of Fig. 
[3]). Viscous disk accretion rate, and hence disk mass, de- 
crease gradually rather than being nearly constant over 
a long time as in the case with a dead zone. If all disks 
have high alpha values a > 10~ 3 , very massive disks may 
be necessary to explain the existence of long-lived, strong 
accretors. 

Fig. [7] shows the same run as Fig. [6j but with the 
planetesimal size of 100 m (RunD2). The planet migrates 
inward, and plunges into the central star within about 2 
Myr. Although the planetary mass increases rapidly as 
the protoplanet migrates, the final mass is about 15 Me, 
and there is little gas accretion. Although the final gas 
accretion rate becomes comparable to the solid accretion 




Fig. 6. — Evolution of planet and disk for planetesimal sizes of 10 
km and no DZ (RunDl). Top: Evolution of a planet with an initial 
mass of 0.6Mb in an evolving disk with no dead zone. The disk 
viscosity parameter is set to a = 2 X 10 — 3 , and the planetesimal size 
is 10 km. We use our fiducial softening parameter B = 0.6. There 
is very little growth of a planet, and as a result, the planet does 
not migrate much. Bottom: Corresponding evolution of planet and 
disk masses as well as disk accretion rate onto the star. Solid, and 
dashed curve shows disk accretion rate, and disk mass, respectively. 
Thick orange, and black curves are on top of each other, and show 
core, and total planetary mass, respectively. 

rate (~5 x 10 -5 Me yr _1 ), migration is too fast in this 
case to make a Jupiter-like planet. It should be noted, 
however, that the resolution of our simulation is good 
down to 0.1 AU, and that we don't define the inner 
dead zone edge. Therefore, any planets which appear 
to be "lost" in our simulations may still be alive if we 
improve our resolution, take account of the inner dead 
zone edge, or include the star-planet tidal interaction. 

Now we show the case of slower migration by adopting 
a larger softening parameter of B = 0.9 (RunD3), which 
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Fig. 7. — Evolution of planet and disk with no DZ (as in Fig. 
6), except with planetesimal sizes of 100 m (RunD2). Top: The 
protoplanet plunges into the star within 2 Myr, and its final mass 
is roughly 20 Me ■ Bottom: Corresponding evolution of planet and 
disk masses as well as the disk accretion rate onto the star. 

roughly reproduces the migration profile of A05 (see their 
§3.1.) Fig. [8] shows such an evolution in a disk with 100 
m-size planetesimals. The planet migrates more slowly, 
but grows only up to ~ IOMe. From disk mass evolu- 
tion (dotted line on the bottom panel), we can see that 
formation of a Jupiter mass planet is impossible after 4 
Myr, since the disk mass falls below Mj. In our run, the 
planet's mass is ~ 7Me at 4 Myr. 

This is very different from what A05 obtained, where 
the planet starts opening a gap after 0.8 Myr at around 
6 AU, and stops its migration at around 5.5 AU. This 
results from the difference in solid accretion prescrip- 
tions as described briefly in the previous section. A05 
followed P96 and adopted rapid planetesimal accretion 
for the growth of the planetary core with an initial mass 




Fig. 8. — Evolution of planet and disk for planetesimal sizes of 
100 m and no DZ, but with larger softening parameter B = 0.9 
(RunD3). The planet migrates much slower, but grows only up to 
9 M E . 

of 0.6Afg. Therefore, their core mass reaches several Me 
by 0.8 Myr, while our core mass is about 2Mb at that 
time. 

In P96, the planetesimal accretion slows when the pro- 
toplanet depletes the planetesimal feeding zone. There- 
fore, during Phase 2, their planetary mass increases at 
the gas accretion rate, which is much smaller than the 
planetesimal accretion rate in Phase 1. As a result, the 
planet formation timescale in P96 was uncomfortably 
long. A05 overcame this problem by constantly replen- 
ishing the planetary feeding zone due to migration, and 
therefore increasing the total mass of a planet at the 
rapid planetesimal accretion rate rather than the gas ac- 
cretion rate. Thus, the planet accretes gas much quicker 
in their case. On the other hand, Fig. [9] shows that the 
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Fig. 9. — Corresponding mass accretion rate for Fig. [8] Orange 
curve is the core accretion rate, and blue curve is the envelope 
accretion rate. As the planetary mass grows, gas accretion rate 
increases and the two accretion rates become comparable to each 
other. 

planetesimal accretion rate in our model is a few orders 
of magnitude lower than that during Phase 1 in P96, 
and comparable to that during Phase 2 in P96. Gas ac- 
cretion rate is even smaller, or at most comparable to 
planetesimal accretion rate throughout the simulation. 
Therefore, in our case, a protoplanet increases its mass 
on nearly constant time scale all the time, independent of 
whether it is migrating or not. In short, for the oligarchic 
core accretion process, we find that the replenishment of 
the planetesimal feeding zone due to migration does not 
help planetary growth as much as it was shown by A05. 

5. PLANET FORMATION IN AN EVOLVING DISK WITH A 
DEAD ZONE 

In this section, we combine all of the elements of proto- 
planetary accretion and disk evolution presented above 
with the dead zone. 

5.1. Without planet migration 

First, we study the case with no migration, but with 
disk evolution. This, of course, is an unlikely scenario, 
because a growing planetary core would quickly migrate 
toward the central star, especially inside the dead zone 
due to an enhanced disk mass (see MPT07). We show 
these cases nevertheless, since they provide reasonable 
estimates of planetary masses achievable in our fiducial 
disks. 

Fig. [TU] and [11] show the results of S = 10 3 and 
10 4 g cm -2 respectively with standard (top panels) and 
reduced (bottom panels) opacities (see §2.3 and Table Q] 
for initial conditions for RunEl(2), and RunFl(2)). The 
core has an initial mass of 0.6Me as before, and stays 
at 5 AU throughout the simulations. The planetesimal 
size is 100 m. In all of these cases, the protoplanetary 
core is originally inside the dead zone, and is left outside 
it as the disk grows and the dead zone shrinks. The fi- 
nal planetary mass of the case with a massive disk and a 
reduced opacity (bottom panel of Fig. ITTj) is about half 
a Jupiter mass, while all the other cases produce plan- 
ets about a Neptune mass (20Me). Note that in both 
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Fig. 10. — Evolution of disk mass and accretion rate, in presence 
of a DZ, and Eo = 10 3 g cm -2 , with standard opacity (top panel, 
RunEl) and reduced opacity (bottom panel, RunE2). Also plot- 
ted are disk mass evolution (dashed curves) along with the disk 
accretion rate onto the central star (thin solid curves.) 

MMSN-mass and more massive disks, the opacity reduc- 
tion leads to more efficient gas accretion as seen in §3.2. 
These results further emphasize the importance of the 
opacity reduction in helping rapid planet formation, and 
indicate that it is difficult to form a Jupiter mass planet 
in a MMSN-mass disk. 

Also plotted in Fig. [TO] and [11] are evolution of disk 
mass (dashed lines) and disk accretion rate onto the cen- 
tral star (thin solid lines.) As shown in §2.2, a dead zone 
enforces a nearly constant accretion rate onto the star, 
and a rapid dispersal of a disk once it's gone. 

5.2. With planet migration 

Now, we finally show the cases with planetary forma- 
tion and migration in an evolving disk with a dead zone. 
The initial conditions are the same as the previous sub- 
section, except that planetary migration is allowed with 
our fiducial softening parameter B = 0.6. 

In Fig. [TOJ and [TOJ we show the evolution of a plane- 
tary core initially at 8, and 10 AU in a MMSN-like disk 
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Fig. 11. — Evolution of disk mass and accretion rate (as in Fig. 
1101 except for column density So = 10 4 g cm -2 , with standard 
(top panel, RunFl) and reduced (bottom panel, RunF2) opacity. 

with S = 10 3 g cm -2 , by assuming a reduced opacity 
and a 100 m-size planetesimal disk (RunGl and G2, re- 
spectively). Since the initial outer dead zone radius is 8.2 
AU, the protoplanet starts inside(outside) the dead zone 
in RunGl(G2). Note that, in Fig. QJ (RunGl), the pro- 
toplanet appears as if it were initially outside the dead 
zone, because the time axis is logarithmic and starts at 
10 5 yr. In this case, the protoplanet is originally located 
inside the dead zone, but is left outside it as the dead 
zone shrinks over time. 

The core building phase takes longer as the initial or- 
bital radius moves outward, and as a result, the outer 
planet accretes less amount of gas. As expected from 
§5.1, neither case leads to a gas giant. Final masses are 
20M E for 8 AU case (RunGl), and 7M E for 10 AU case 
(RunG2). In RunG2, the core accretion is so slow that 
the planet barely migrates compared to RunGl. We ob- 
tain a similar result to RunGl for an initial radius just 
outside the dead zone at 9 AU. 

From the comparison of Fig. [T2]with Fig. the effect 
of a dead zone is apparent. Instead of losing a growing 
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Fig. 12. — Evolution of planet and disk with a planetary core 
initially just inside DZ at 8 AU. The disk column density So = 
10 3 g cm -2 , reduced opacity, and a 100 m-sizc planetesimal disk 
are assumed (RunGl). Top: Evolution of both disk and planetary 
orbital radius. Note that the planet appears as if it were initially 
outside the DZ, because the time axis is logarithmic and starts at 
10 5 yr. Bottom: Corresponding evolution of planetary and disk 
masses as well as the disk accretion rate onto the central star. 

protoplanet to the central star, the dead zone stops the 
migration of the protoplanet by balancing the inner and 
outer torques. As a result, the protoplanet has a chance 
of growing further at a slower, viscous accretion rate of 
the disk. 

Fig. HH and [15] show the evolution of a planetary 
core initially at 15 AU in a more massive disk with 
Eo = 10 4 g cm~ 2 , by assuming a reduced (RunHl), and a 
standard (RunH2) opacity, respectively. For both cases, 
a 100 m-size planetesimal disk is used. Since the initial 
outer dead zone radius is 15.8 AU, these planets start 
just inside the dead zone. However, the protoplanets 
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Fig. 13. — Evolution of planet and disk as in Fig. 1121 except the 
initial orbital radius of a planet is outside DZ at 10 AU (RunG2). 

appear as if they were initially outside the dead zone, 
because the time axis is logarithmic and starts at 10 5 yr. 
Run HI obtains a gas giant with ~ 0.3Afj in roughly 2.5 
Myr, while in Run H2, the protoplanet grows only up 
to ~ 20Mb. This demonstrates that we need an opacity 
reduction (or some other mechanism) even when we in- 
clude planet migration. As emphasized in §4, this is the 
result of a slower oligarchic core accretion compared to a 
runaway core accretion. The overall mass accretion onto 
a protoplanet does not speed up as a result of migration, 
and hence the constant replenishment of the planetesimal 
feeding zone. 

Also, by comparing the bottom panels of Fig. [T3] and 
1151 it appears that a gap-opening planet speeds up the 
dispersal of a gas disk. This additional disk dispersal 
mechanism should be investigated in multiple planet for- 
mation models. 
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Time [years] 

Fig. 14. — Evolution of planet and disk as in Fig. 1121 except 
disk column density of So = 10 4 g cm -2 , and the initial orbital 
radius of a planet just inside DZ at 15 AU (RunHl). Note that the 
protoplanet appears as if it were initially outside the DZ, because 
the time axis is logarithmic and starts at 10 yr. The second jump 
in mass seen on the bottom panel occurs at semi-major axis smaller 
than our resolution limit, and may not represent the true evolution 
of the planet. 

6. DISCUSSION AND CONCLUSIONS 

We have presented some new results on planet forma- 
tion and migration in evolving disks with dead zones. 
The most significant is that dead zones provide a natural 
way of saving planetary systems even as the planets mi- 
grate through disks whose properties change significantly 
over hundreds of thousands to millions of years. The 
dissipation of the disk does place interesting constraints 
on their masses - we find that only Neptunian mass 
planets can be formed in the MMSN-mass disk mod- 
els. Jovian planet formation requires more massive disks, 
which has also been suggested by other groups (e.g. P96 
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Fig. 15. — Evolution of planet and disk as in Fig. [H but the 
standard opacity instead of reduced opacity is assumed (RunH2). 
Again, although the initial protoplanetary radius is just inside DZ 
at 15 AU, the protoplanet appears as if it were initially outside the 
DZ, because the time axis is logarithmic and starts at 10 yr. 

and HBL05), and was recently d emonstrated by multipl e 
planet formation simulations by iThommes et al.l (|2008f) . 
The time scale that we find for Jovian planet formation 
in these more massive disks - about 2.5 Myr - is within 
observational limits of disk lifetimes. While not drasti- 
cally reduced from previous estimates, our Jovian planet 
formation time incorporates many new aspects including 
migration in the presence of dead zones, effects of slower 
oligarchic growth that are more realistic than earlier ac- 
cretion scenarios, and an effect of a reduced opacity. 

Dead zones turn out to play a potentially important 
role as blockades against inward planetary motion. Our 
planetary cores are generally kept outside of the dead 
zone - and therefore immersed in the region where they 
can "feed" more effectively on surrounding gas. When 



the planets become massive enough to open a gap, the 
outwards directed torque associated with the density gra- 
dient of the outer edge of the dead zone weakens, and 
the planets migrate into the dead zone. Planets migrate 
slowly inside the dead zone due to the lower viscosity 
there compared to outside it. 

The methods we used were straightforward to imple- 
ment in our planetary/disk evolution code. In order to 
model dead zone evolution, we simply assumed that the 
dead zone is where the surface mass density is above 
the critical value (S > E cr i f = 21, and 80 g cm~ 2 for 
So = 10 3 and 10 4 g cm -2 respectively) that was de- 
termined from a stationary disk model in MP06. Also, 
we have included gas accretion through the surface lay- 
ers, which contributes to make a dead zone shrink very 
rapidly (see Fig. [3|). Combining these two conditions, 
we have found that the dead zone radius shrinks from 
8.2 AU to 2 AU within 2 Myr for E = 10 3 g cm" 2 . 

For the planet formation part of the code, we fol- 
low the approach by P96, but with a few differences. 
First of all, we assume that the planetesimal accretion 
stage is mainly carried out by oligarchic growth, rather 
than by rapid runaway accretion used by P96. This is 
a reasonable assumption since the runaway growth most 
likely ceases long before the planetary m ass reaches a few 
times 10~ 5 Me (jThommes et al.ll2003l ). which is much 
smaller than our initial planetary mass (0.6Mf;), and 
switches to a slower oliga rchic growth (|Ida fc Makinol 
Il993t iKokubo k Idalfl998f l. Depending on the size dis- 
tribution of planetesimals, the rate of accretion could 
be significantly different. Secondly, we p arameterize our 
;as ac cretion rate following INEOO and iDAngelo et al.l 
2003). Therefore, we are not calculating planetesimal 
and gas accretion rates in an interactive way as in P96, 
HBL05, or A05. However, our simulations show reason- 
able agreements with their results (see §3). Also, since 
we take account of the subdisk accretion phase, the final 
stage of gas accretion slows down, rather than exponen- 
tially increases as in P96, or artificially cuts off as in 
HBL05. Thirdly, we include a ID gas disk evolution and 
planet migration. Our approach is similar to A05, but 
we don't include photoevaporation effects for this study, 
since this is likely to be negligible during planet forma- 
tion. 

There are several aspects of planet formation that we 
have not included explicitly in this study. We list several 
below, but note that we do not expect their absence to 
strongly affect our results. 

First of all, our simple torque prescription does not 
capture the nature of a turbulent disk properly. A 
number of numerical simulations have shown that tur- 
bulent fluctuations c an cause torque fluctuatio ns (e.g. 
lLaughlin efail 12004 iNelson fc Papaloizoull2004h. which 
leads to stochastic type I migration (jJohnson et alJ 
2006). This results in even faster migration for most 
planet s, but also allows o utward migration for some of 
them (| Johnson et al.ll2~006[ K The inclusion of such an ef- 
fect may allow a planet outside the dead zone to "jump 
the barrier" provided by the edge of a dead zone. Even 
within a dead zone, turbulent fluctuations in upper and 
lower active layers may be able to affect planet migration 
significantly. 

Also, we do not take account of the evolution of tem- 
perature profile of a gas disk. To better evaluate the disk 
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evolution, we should include the radiative transfer to the 
code. This would affect not only the dead zone and disk 
evolution, but also planet migra tion. A recent work of 
iPaardekooper fc Mellema (2008) demonstrated that, in 
a non-isothermal disk, planet migration is preferentially 
outward at around 5 AU until the disk mass decreases 
significantly. Therefore, more precise treatment of a disk 
temperature may save type-I migrators effectively, even 
inside a dead zone. 

Perhaps the most important simplification of our 
model is the ID treatment of the disk, and the rather 
crude prescriptions for the dead zone, which led to a 
very sharp density transition at the outer dead zone ra- 
dius (see Fig. [3]). This feature enabled a planet to stop 
its migration by balancing inner and outer torques. How- 
ever, such a sharp density gradient makes a disk locally 
Rayleigh unstable (i.e. an epicyclic frequency n 2 < 0), 
which may render the sharp transition from the dead 
zone to the outer active disk much smoother. Even when 
the disk is still locally stable, the Rossby wave instabil- 
ity can smooth the density gradient in a similar manner 
when pressure varies signific antly over a few times the 
disk thickness (|Li et al.ll2001l . and the references therein) . 
This possibility should be checked carefully by numerical 
simulations. 

We did not take account of the effect of the corota- 
tion torque either. This is a safe assumption most of the 
time, since the corotation torque depends on the gradi- 
ent of the inverse of the specific vorticity (<9(E/ B) /dr) 
(|Goldreich fc Tremaind Il980b iMassetl l200lh . and there- 
fore its effect becomes particularly important when the 
surface mass density changes sharply, e.g. at the edge 
of a dead zone. However, we do not expect this affects 
our results significantly, because inner and outer Lind- 
blad torques balance with each other far away from the 
density jump, where the corotation torque is likely neg- 
ligible, and planet migration stalls (also see Appendix B 
in MPT07). When Rayleigh, or Rossby wave instability 
makes the density gradient at the outer edge of a dead 
zone smoother, the corotation torque can become impor- 
tant. In such planet would not be stopped at 
the outer edge of a dead zone, but pulled into it instead. 

The corotation torque is also found to be important 
for migration of sub- Jovian mass planets like Saturn 
(|Masset fe Papaloizoul [20031 . In a standard disk, the 
time scale of the so-called type III migration for such 
planets is comparable to, or even sho rter than, type I 
migration (M asset fc Papalo izou 2003). Although this is 
most relevant to our RunHl (Fig. [14]), we don't expect a 
significant change in our result. This is because type III 
migration is likely to switch to slower type II migration 
as soon as the planet enters a dead zone, where planets 
tend to open a wider gap at smaller mass due to its low 
viscosity. 

Another simplifications is that we have only consid- 
ered a protoplanetary disk with a single core, while the 
real disks are expected to have mul tiple cores . Stud - 
ies of this kind have been done by IChambersI ([2006) ; 
iThommes et al.l (I2008T) for disks with no dead zone, and 
by iMorbidelli et alj (|2008f) at the in ner edges o f disk s 
with dead zones. A recent paper by llda fc Linl (2008), 
took a similar approach to us, and studied a dead zone's 
effects on retaining icy grains, as well as protoplane- 
tary cores. They reproduced the observed frequency 



and mass-period distribution of gas giants around solar- 
type stars with a moderate reduction in type I migration 
speed. 

Yet another simplification is that all of our planets in 
this study are on circular orbits. Planet-disk interaction 
has been proposed as a method to drive planetary ec- 
centricity (|Goldreich fc Sarill2003[ ), and recent numerical 
simu lations confirmed this for massive (> M A pla nets 
(e.g. iMasset fc Qgjfviel l200i iD'Angelo etall 12001 . If 
a planetary eccentricity is enhanced significantly (e 3> 
h /r) , then planet migration could be slowed down 
(e.g. iPapaloizou fc Larwoodl[2000l : iPapaloizou 2002 ), and 
mass growth rate could in crease (e.g. ID'Angelo et al.l 
[200llKlev fc Dirksenl[200l . 

We also employ simplified planetesimal disks. First, we 
only consider single-size planetesimals ranging from 100 
m to 100 km, while the real disks should have multiple- 
size planetesimals. Secondly, we don't include the plan- 
etesimals' effect on pl anet migration, whic h could poten- 
tially be important ([Murray et al.l ll998). Thirdly, we 
assume the dispersion-dominated random velocities for 
planetesimals with all sizes. Since small planetesimals 
experience strong gas drag and achieve reduced random 
velocities, their mass accretio n rate should be treated 
in th e shear-dominated regime (|Rafikovil200l iChamberd 
[2Q0l). This may be im portant for 100 m-size planetes- 
imals ([Chambers! [2006) . and the planetesimal accretion 
would be runaway, rather than oligarchic in such a case. 

Finally , we don't include the effects of photoe vapora- 
tion (e.g. lShu et aTlflliM iHollenbach et al.lll994h . which 
is likel y to be important during the last stage of d isk evo- 
lution (|Clarke et al.ll2001l ; [Alexander et al.ll2006f) . How- 
ever, photoevaporation becomes important only when 
disk mass becomes significantly low (~ lMj or so), and 
therefore it is unlikely to affect our results. 

We conclude by listing our major findings: 

(1) Dead zones strongly evolve over the duration of 
the disk, starting with outer radii of order 10 — 15 AU, 
and shrinking with time to of order an AU or less. 

(2) Protoplanets which are left outside the dead zone 
tend to migrate inward as they accrete planetesimals, 
and stop just outside the outer dead zone radius due to a 
steep surface mass density jump. Then such protoplan- 
ets migrate at the viscous time scale of the shrinking 
dead zone until they achieve gap-opening mass. Once 
they open a gap, they enter the dead zone, and mi- 
grate inwards rather slowly due to the low viscosity there. 

(3) The final mass of a planet is determined by disk 
mass. In a minimum-mass solar nebula disk, it is 
difficult to get a planet more massive than Neptune. 
We find Jovian planets form in disks that are ten times 
more massive than this. 

(4) A Jupiter mass planet can form within ~ 2.5 
Myr (see Fig. [T4|) in a disk with a dead zone, by 
assuming standard type I migration. However, since 
planet migration does not help core accretion as much 
in oligarchic growth as in runaway growth scenario, 
we expect that the opacity reduction (or some other 
mechanism) is necessary to form a Jovian planet within 
a disk life time. 
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(5) Dead zones may help explain the existence of 
long-lived strong accretors. When a dead zone is 
present, the mass accretion rate is likely to take a nearly 
constant value until the dead zone disappears, rather 
than decreasing gradually as expected in a standard 
disk without a dead zone. Once this happens, the rest 
of the disk is dispersed on a few Myr time scale, unless 
photoevaporation and/or planet formation provide 
additional sources of disk dissipation. 
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